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A large set of S'pitzer-selected transitional disks in the Ophiuchus molecular cloud was recently examined by Cieza et 
■ al (2010), and 4 of the targets were identified as (giant) planet-forming candidates based on the morphology of their 

[•-p'] ' Spectral Energy Distributions (SEDs), the apparent lack of stellar companions, and evidence for accretion. Here we 

\ characterize the structures of these disks modeling their optical, infrared and (sub)millimeter SEDs. We use the Monte 

, Carlo radiative transfer package RADMC to construct a parametric model of the dust distribution in a flared disk with 

an inner cavity and calculate the temperature structure consistent with the density profile, in thermal equilibrium with 
the irradiating star. For each object, we conducted a Bayesian exploration of the parameter space generating Monte 
Carlo Markov Chains (MCMC) that allow the identification of the best fit parameters and to constrain their range of 
statistical confidence. Our calculations point to the presence of evacuated cavities with radii ~ 2 — 8 AU, consistent with 
^ ■ having been carved by embedded giant planets. We found parameter values consistent with those previously given in 

^ ' the literature, indicating a mild degree of grain growth and dust settling, which deserves to be investigated with further 

modeling and follow up observations. Resolved images with (sub) millimeter interferometers are required to break some 
of the degeneracies of the models and better constrain the physical properties of these fascinating disks. 

> 

Xf^ 1. Introduction Probing the structure of disks that are suspected to 

' be forming planets is the most promising approach to 

I The infrared spectral energy distributions (SEDs) of understand the conditions in which planets are formed. 

circumstellar transitional disks reveal the presence of an The best indication for ongoing planet formation in disks 
optically thin inner region and an optically thick outer is the detection of tidal gaps (e.g. Pietu et al. 2006) 
disk. Several mechanisms relevant to the overall evolution corresponding to a ring with significant decrease in the 
_ ^ ^ of circumstellar disks, and in particular to the short-lived surface density (or the whole inner disk if it was depleted 
\ phase when they dissipate, have been proposed to explain by accretion, e.g. Varneiere et al 2006). A spectacular 
. the so-called opacity holes of transition disks: giant planet confirmation has arrived with the recent detection of 
C3 \ formation, grain growth, photoevaporation, and tidal the first potential substellar object within the gap of 
■ ■ ■ ' truncation in close binaries. See Williams k Cieza, 2011 the transitional disk T Chamaleontis (Huelamo et al. 
for a recent review. The processes responsible for the inner 2011). Inner holes and gaps have already been observed 
holes of transition disks can tentatively be distinguished at (sub) millimeter wavelengths in a handful of objects 
when disk masses, accretion rates, and muhiplicity in- bright enough for resolving disk structure (Pietu et al. 
formation are available (Najita et al. 2007; Cieza 2008). 2006; Hughes et al. 2007, 2009; Brown et al. 2008, 2009; 
Following this approach, Cieza et al. (2010, hereafter Andrews et al. 2009, 2010a, 2011, Isella et al 2010a, 2010b). 
Paper I) presented the initial results of an ongoing project 
to characterize a large set of transition disk candidates in 

nearby star-forming regions. ^ , . , c ■ . ■,• j- i 

o o Comprehensive studies of similar transition disks are 

necessary to increase the empirical constraints on their 

structures and investigate their diversity. This is the pri- 

Send offprint requests to: Mariana Orellana morel- mary motivation of this work, where we apply a paramet- 

lana@fcaglp.unlp.edu.ar ric description to the 4 planet- forming disks candidates in 
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Ophiuchus presented in Paper I and derive the best fitting 
physical characteristics for their known SEDs. 

2. The targets 

In Paper I, the observed SEDs were characterized by two 
parameters, as introduced by Cieza et al. (2007): the longest 
wavelength at which the observed flux is dominated by the 
stellar photosphere, Aturn-off, and the slope of the IR ex- 
cess, aoxccss, computed between Aturn-off and 24 jim. The 
former parameter correlates with the size of the inner hole, 
and the latter with the sharpness of the edge of the hole, 
i.e., a large increase in the dust density over a small range 
in radii is indicated by a positive aexcess- Accreting disks 
with sharp inner holes that seem to lack stellar compan- 
ions (e.g., from adaptive optics imaging) stand as the most 
promising candidates for ongoing (giant) planet formation. 

There were 4 cases fulfilling these criteria in Paper I 
from a sample of 26 of Spitzer-selected transitional disks 
in the Ophiuchus molecular cloud {d ~ 125 pc, Loinard 
et al. 2008). Their basic properties are listed in Tabled] 
The mass accretion rates for our targets, inferred from the 
Ha emission line widths (Natta et al. 2004), range from 
1.3x10-1° to 5x10-8 Mo/yr. 

For all of them we have at hand the photometric fluxes 
in the i?-band from the USNO-Bl catalog, the J, i?, and 
Ks bands, from the Two Micron All Sky Survey (2MASS), 
as well as fluxes at 3.6, 4.5, 5.8, 8.0, 24, and 70 /im obtained 
by Spitzei^. We have assumed errors that are within obser- 
vational standards (^ 20%) . In addition, all of our targets 
have been observed with the Submillimeter Array (SMA). 
For two of them (Tran 21 and 32) we obtained the fluxes 
at 1.3 mm, while for Tran 11 and 31 we could only de- 
rive upper limits (see Paper I for details on SED data). For 
Tran 32, Fgso^m was given in Nutter et al. (2006). One of 
our targets, Tran 11, was recently included in the study 
of cold disks (disks with large inner dust holes) performed 
by Merm et al (2010), i.e. their source ^ 24. It was mod- 
eled including a positive detection from MIPS photometry: 
^70Aim — 537 ± 78.8 mJy, and the measured Spitzer-IRS 
spectrum, but no constraints of the flux at (sub)-millimeter 
wavelengths. For the sake of completeness, we have redone 
the fit to the SED of this source with our set of free param- 
eters and including the 1.3 mm fiux upper limit. 

3. Radiation transfer models 

The modeling that we have applied is similar to those per- 
formed by Andrews et al. (2009, 2010a) and Brown et 
al. (2009), who have confirmed their physical estimates 
from SED modeling through direct imaging. A 2-D struc- 
ture model for fiared disks is combined with the Monte 
Carlo continuum radiative transfer package RADMC v3.1 
(Dullemond & Dominik, 2004), modified to include a den- 
sity reduction as an inner cavity. The code computes a tem- 
perature structure consistent with the given density profile, 
and in equilibrium with the irradiation by the central star. 
The disk is presumed to be passive, an assumption that is 
supported by the low disk to stellar luminosity ratios of our 
sample, 550.005 according to estimates in Paper I. 
We consider a surface density profile characterized by a 

^ two of the 70 fim measurements given in Paper I are upper 
limits. 



power-law, S oc i? with an exponential taper at larger 
radii (cx e"^^/^''^ ^), where Rc is the characteristic radius. 
This is physically motivated by the success of similarity so- 
lutions of viscous disks to reproduce the observed gradual 
density decay at large radii (Hughes et al. 2008). S is nor- 
malized to obtain the total mass of the disk, Md, when 
integrated. The radial index was fixed to be 7 = 1 which is 
a typical value within the range 7 = 0.4 — 1.1 established 
by Andrews et al (2010a). Our option could be questioned 
in the light of results by Isella et al. (2009) who have in- 
dependently inferred slopes from steep to quite shallow in 
their sample of spatially resolved disks. 
Resolved images are therefore mandatory to obtain more 
accurate estimates of 7 for our particular targets. We set 
the characteristic radius Rc = 100 AU. However, there is 
no spatially resolved information in the SEDs alone, and 
the data can be reproduced equally well with a wide range 
of outer disk values (Andrews et al. 2010b). The value we 
choose is representative of the disks with resolved interfer- 
ometric visibilities, that are Rc = li — 198 AU (Andrews et 
al. 2009, 2010b) in Ophiuchus, and i?c ^ 30 - 230 AU for 
Taurus- Auriga (Isella, Carpenter & Sargent, 2009). Note 
that larger outer radii (100 - 1100 AU) have been obtained 
with different fitting techniques (sharply truncated power 
law fits to CO observations) and are not directly compara- 
ble with the Rc values (see Williams & Cieza 2011). Aside 
from the extreme case of a nearly edge-on viewing angle, 
the disk inclination cannot be determined from unresolved 
observations. Scattered light images have proben useful in 
this sense (Pinte et al. 2007). We have set an intermediate 
representative inclination i — 30° in our modeling. 
There is strong observational evidence that circumstellar 
dust can present some degree of settling to the midplane 
(Furlan et al. 2006, McClure et al. 2010), which is in agree- 
ment with theoretical predictions (e.g., Dominik et al. 2007 
and references). Indeed, the growth of the dust grains, 
which proceeds in a complex fashion, can be accelerated 
by dust settling, and it is expected that the larger particles 
are aggregated close to the equatorial plane. The dust can 
therefore be distributed differently than the gas, which is 
expected to remain in hydrostatic equilibrium (Chiang & 
Goldreich, 1997). This is described by a vertical Gaussian 
density profile with scale-height H and radial index ip, i.e. 
H (X R^+'l', anchored at i? = 100 AU by iJioo- Vertically 
extended disks (with larger H) are heated more efficiently 
and thus reemit more radiation at mid- and far-IR wave- 
lengths. Here we have fixed ip = 0.2 in order to reduce the 
number of free parameters of the models. 
Following Andrews et al. (2009), we consider a radius, i?cav, 
such that the modified surface density is E' = (5cavS when 
R < -Rcav, and Scav < 1 artificially reduces the density in- 
side the cavity, in order to mimic the deficit in the disk 
emission. The whole disk extends down to the distance. 
Rim where the temperature is enough to sublimate the 
dust grains (~ 1500 K). Our 4 targets present a small or 
negligible NIR excess above the photospheric value (note 
the Aturn-off values in Table [T]), suggesting cavities that are 
largely evacuated of dust, therefore we have varied log((5cav) 
as a sensitive parameter to be established by the data. 
We have adopted the opacity spectrum given by Andrews 
et al. (2009) who considered the silicate and graphite abun- 
dances determined for the ISM, and updated optical prop- 
erties. The grain size distribution, n{a) oc a^^'^, extends in 
diameters from a^i^ = 0.005 /xm up to amax=l mm, and 
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Table 1. Target Properties 



Name" 


Spitzer ID 


SpT 










Age 


Av 


log Mfcc 


Aturn — off 


Q^cxccss 






SSTc2d_ 




(K) 


(Le) 


(Re) 




(Myr) 


(mag) 


(Me/yr) 


(Mm) 




(Mjup) 


Tran 11 


J162506.9-235050 


M3 


3470 


0.24 


1.25 


0.3 


2.1 


3.8 


-8.8 


8.0 


0.65 


<1.5 


Tran 21 


J162854. 1-244744 


M2 


3580 


0.51 


1.74 


0.4 


1.4 


5.4 


-9.3 


8.0 


0.69 


1.3 


Tran 31 


J163205.5-250236 


M2 


3580 


0.19 


1.08 


0.4 


4.1 


5.0 


-7.3 


8.0 


0.30 


<1.3 


Tran 32 


J163355. 6- 244205 


K7 


4000 


0.78 


1.70 


0.7 


2.0 


5.0 


-9.9 


8.0 


0.72 


11.1 



" Alternative names: Tran 21 is WSB 63, Tran 31 is WSB 75, and Tran 32 is RXJ1633. 9-2242. Stellar parameters from pre-main 
sequence evolutionary tracks by Siess et al (2000). Based on the velocity dispersion of the Ha line from Paper I. Rough estimates 
of the disk masses based on a single (sub)millimeter flux or upper limits from Paper I. 



the opacities were calculated with a simple code for Mie 
scattering. 

The stellar properties given in Table [1] are in each case 
fixed inputs for the models, and Kurucz spectra are used 
as models of the central stars. As the disk SED is not in- 
fluenced by the A 5s 3 /im fluxes, they were not included in 
the final fit, avoiding a systematic error term (shift in the 
X^)- For each disk, we have explored the parameter space, 
{Mdisk, Hioo, i?cav, log Jcav}, by generation of Monte Carlo 
Markov Chains, which is a parameter space exploration 
technique designed to provide the best fitting values and 
their uncertainties. This technique attains a better perfor- 
mance than classical data fitting methods (e.g., Press et al. 
1992; Gregory 2005) for problems with large numbers of pa- 
rameters, which is our case. In the fitting procedure we have 
followed Ford (2005) and Gregory (2005). The marginal 
probability distribution for each parameter where derived 
from the MCMC chains. The 90% percentile of such distri- 
butionf0 simmetrically centered in the mode when possible 
(or using the larger side if not) was used to define the error 
estimates within a conservative approach. 



to flat geometries (-ffioo ^ 4 AU). However, the on-line fit- 
ting tool fails to accurately reproduce the SED for 2 of our 
targets, Tran 11 and 32. A discussion on the preferential 
use of RADMG to model transitional disks and comments 
on alternative observational characterization of them can 
be found in Merfn et al (2010). 



Table 2. Model parameters. 



# 



(Mjup) 



-"100 

(AU) 



RcELV 

(AU) 



log Sea. 

(dex) 



11 


<2.0 ± 0.7 


<4.2 ± 1.3 


4.8 ± 2.5 


^ -6.2 


21 


0.6 ± 0.2 


3.1 ± 0.8 


1.9 ± 0.3 


-4.9 ± 0.8 


31 


< 1.7 ± 0.1 


<1.9 ± 0.1 


1.5 ± 0.4 


-5.2 ± 0.5 


32 


17.0 ± 5.3 


2.0 ± 0.5 


7.9 ± 2.3 


^ -6.3 



" Mdisk are upper limits for Tran 11 and 31 as their 1.3 mm 
fluxes are also upper limits. '' Hioo is an upper limit for Tran 31 
as its 70 fim fluxes is also an upper limit. 



4. Results 

The estimated parameters that best reproduce the data for 
our planet forming disk candidates are listed in Table [21 
while the corresponding SEDs are shown in Figure [1] These 
fits have reduced — 0.9 — 10, i.e. the considering 
only the observational data point errors, and divided by 
the number of degrees of freedom. The disk masses ob- 
tained are within a factor ~2 than the rough estimates from 
Paper I, which were based on a single (sub)millimeter flux. 
We find iJioo ^2-6 AU in all four cases, with estimated 
uncertainties i.e. ^ 1.5 AU. The comparison between the 
dust scale heights to the scale heights of the gas, which is 
in hydrostatic equilibrium, give a settling ratio 0.13-0.25. 
Specifically for Tran 32 we get a rather flat geometry, with 
the smallest i?ioo in our sample. This disk presents also the 
largest cavity and is bright enough to be detected by SMA 
extended configuration. A simultaneous SED -I- image fit 
will be presented in a future paper. 

To test our results against an alternative model, we 
have used the precomputed grid of SEDs by Robitaille et 
al. (2007). A cavity is also inferred in all 4 cases, with 
.Rcav ~ 4 — 38 AU, as well as some dust settling leading 

^ Note that there is no guarantee that the marginal distribu- 
tion is Gaussian so the dispersion of the sample can only be a 
rough indicator of the probable error as in Isella et al (2009). 



5. Discussion 

Extensive work in non-linear model fitting have been ap- 
plied in many areas, and recently to disk modeling by Isella 
et al. (2009). They have used the two layer approxima- 
tion of Chiang & Goldreich (1997), whereas we compute 
the thermal structures with significant detail using the 
RADMC code. In this study we have followed the litera- 
ture by Andrews et al. (same parametric model and opac- 
ities) who have performed refined grid searches over 8 pa- 
rameters. Their acute comments on the technical obstacles 
(Andrews et al. 2011) and non-imiqueness of the fit, also 
apply for our approach. 

We have estimated some of the physical parameters de- 
scribing the dust structures of 4 transition disk systems that 
are excellent candidates for ongoing giant planet formation. 
The cavities sizes inferred for our targets, iicav = 2 — 8 AU, 
are in agreement with the distribution of semi-major axis 
of exoplanets, which has a bi-modal behavior, peaking at 
~ 0.05 and 2 AU and extending up to ~ 10 AU (exoplan- 
ets database, |http://exopla nets. org). The mass of the disks 
and the actual size of their inner holes depend on the as- 
sumed opacities. The small scale height of the dust in all 
our targets suggest this could be a characteristic property 
of planet forming disks. This is an intriguing results that 
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Fig. 1. Best fit SEDs (solid lines) obtained using RADMC, 
in the context of a flared disk model with an inner cavity. 
The dashed line shows the input stellar model atmosphere. 
The filled circles are extinction corrected values, while the 
arrows are upper limits. 



deserved to be investigated with further modeling and fol- 
low up observations. 

An effort was done here 10'' model runs for disk) to 
provide the set of best fitting parameters and their confi- 
dence regions through the application of Bayesian methods, 
but the real weight of the results could be yet unclear given 
the degeneracy of the model into some of the parameters we 
have fixed. An illustration is the case of Tran 11, for which 
a cavity radius i?cav = 3 ± 2 AU, and a settling parameter 
that translates into -0 = 0.07 have been found by Merfn et 
al. (2010) when setting Rc = 200 AU. We used = 100 
AU, V = 0.2 and obtained i?cav = 4.8 ± 2.5 AU. 

The disk masses obtained here cover a rather wide range 
of masses, from ~ O.GMjup to ^ 17Afjup that extends to 
lower values the mass range obtained by Andrews et al., 
i.e. 5 — 42Mjup for systems that have been modeled in 
the same way, i.e. with same parameterizations and opacity 
tables. However, Andrews et al. (2009, 2010a) fitted their 
systems with larger values of the depletion factor, (Scav ~ 
10"''- 10"^, and obtain larger cavities, i?cav = 20-40AU. 
We note that in Andrews et al. (2011) where the model 
include a more complex surface density profile (i.e. with an 
inner disk inside the cavity) the range of masses estimated 
for 12 disks with resolved images is ^ 8 — 128Afjup but com- 
parisons could be misleading in this case. Some differences 



are probably the result of a selection effect. While Andrews 
et al. selected their targets based on large (sub)millimeter 
fluxes, our selection criterion is based on the slope of IR 
SEDs, and the presence of accretion. 

At a southern declination of '^25 deg and a distance 
of ~125 pc, our targets are excellent targets for follow up 
studies with Herschel and the Atacama Large Millimeter 
Array (ALMA) to investigate their properties in more de- 
tail, and in particular, the small holed transition disks will 
require its better spatial resolution. Herschel far-IR pho- 
tometry would bridge the gap between the mid-IR and the 
(sub)mm wavelengths accessible from the ground and help 
to better constrain the scale heights and the flaring angles 
of the disks. Herschel spectroscopy of fine structure lines, 
such as the 63.2 /xm [O I] line, could help to probe their gas 
content and hence the gas to dust mass ratio. Similarly, 
resolved images with ALMA will break some of the degen- 
eracies of the models and will allow to better understand 
the physical properties of these fascinating disks, thereby 
helping to elucidate the conditions in which planets are 
formed. 
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